Theory of non- Fermi liquid near a diagonal electronic nematic state on a square lattice 
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We study effects of Fermi surface fluctuations on a single-particle life time near the diagonal electronic 
nematic phase on a two-dimensional square lattice. It has been shown that there exists a quantum critical point 
(QCP) between the diagonal nematic and isotropic phases. 1 We study the longitudinal fluctuations of the order 
parameter near the critical point, where the singular forward scattering leads to a non-Fermi liquid behavior 
over the whole Fermi surface except along the k x - and fc^-directions. We will also discuss the temperature and 
chemical potential dependence of the single-particle decay rate. 



Introduction — In strongly correlated electronic materials, 
electrons can organize themselves into a pattern with period- 
icity differing from that of the underlying lattice along a par- 
ticular direction. The phase with such a non-trivial charge 
ordering breaks a discrete translational symmetry of underly- 
ing lattice along one direction, and has been named as stripe. 
The stripe physics has been of great interest to both theo- 
retical and experimental condensed matter physicists, espe- 
cially because its quantum fluctuation may play an impor- 
tant role in understanding the mechanism of high tempera- 
ture superconductivity. 2 ' 3,4 The clear evidence of such a novel 
paired stripe order was reported by Mori et al in electron 
diffraction experiment on (La,Ca)Mn03. 5 

It has been proposed that phases with fluctuating stripes can 
be generic ground states of a doped Mott insulator introduced 
by the hole doping. 6 ' 7 These intermediate forms of matter have 
been dubbed as electronic smectic and nematic phases. While 
the smectic breaks a translational symmetry, the electronic ne- 
matic phase breaks a discrete rotational symmetry of underly- 
ing lattice, which can be viewed as fluctuating stripes. The 
electronic nematic phase has been referred as Pomeranchuk 
instability of Fermi liquid, 8 and has been studied in theoreti- 
cal models such as t-J, (extended) Hubbard, and Fermi liquid 
models i 9 ' 10 i n ' 12 The electronic nematic phase has been sug- 
gested as a phase responsible for novel phenomena observed 
in strongly correlated systems. These phenomena include, for 
example, the anisotropic resistivities in high Landau levels of 
quantum Hall systems^ the anisotropic scattering patterns 
observed in one of high temperature superconductors,-^^ 
and the metamagnetic transition and anisotropic transports in 
Sr 3 Ru 2 7 «i£ 

In addition to the anisotropic nature of electronic nematic 
phase, it was shown that the gapless Goldstone mode of the 
nematic phase in continum model leads to the non-Fermi liq- 
uid behavior in single-particle scattering rate on Fermi sur- 
face except along the symmetric directions. 19 However, us- 
ing the same model of quadrupole density interaction on the 
square lattice, it was shown that the putative continuous tran- 
sition is preempted by a first order transition due to the van 
Hove singularity. The regime of nematic phase diminishes 
exponentially as one approaches the continuum limit of small 
density with approximately circular Fermi surface. 20j — There- 
fore, there is no critical Fermi surface fluctuations for the 



non-Fermi liquid behavior in realistic systems with underly- 
ing lattice; both transverse (phase) and longitudinal (ampli- 
tude) modes are gapped. On the other hand, it was later found 
that additional uniform interaction in the forward scattering 
channel suppresses the first order transition. 22 This restores 
the non-Fermi liquid behavior at the Fermi surface, except 
along the zone diagonal direction, via the coupling to the fluc- 
tuation of order parameter at the quantum critical point. 23 

In this paper, we present theory of non-Fermi liquid near 
the diagonal nematic phase. It is worthwhile to note that there 
are two distinct nematic states in the square lattice. 1 Previous 
studies provide understanding of phenomena related to Fermi 
surface distortion along the parallel direction of k x - or k y - 
axis of the Brilluoin zone. The diagonal nematic phase can be 
characterized by a Fermi surface distortion along the diagonal 
direction of the Brilluoin zone. It is shown that the transition 
between the diagonal nematic and isotropic phases is second 
order, because it does not involve the splitting of the van Hove 
singularity in the density of states,' We find the non-Fermi liq- 
uid behavior in the single-particle life time at the Fermi sur- 
face , through the coupling to the longitudinal mode of diag- 
onal nematic order, except along the zone parallel directions 
(x- and y-axis of Brillioun zone). Our numerical results of self 
energy show the expected exponent of 2/3 in frequency simi- 
lar to the parallel nematic phase 23 and the continuum model 19 . 
However, the physical origin of the non-Fermi liquid behavior 
is different for these two cases. 

Model — The model Hamiltonian for the parallel and the 
diagonal nematic phases on a square lattice is written as fol- 
lows, 
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where Ci(k) = cosfc x — cosk y , and C 2 (k) = 2 sin k x sin k y 
are the form factors originating from the quadrupole density- 
density interaction 19 , and £k is the tight binding dispersion 
£k = — 2t(cosk x + cosky) — \L. We assume short range 
quadrupolar interactions of F 2 (q) and G 2 (q) with the form 
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FIG. 1: The spectral function of the effective interaction, Imx/v, as 
a function of uj at slightly away from the QCP for several choices of 
q along (a) zone parallel (0, tt) and (b) zone diagonal (it, tt) direc- 
tions. There is a sharp propagating mode along the (0, n) direction, 
and the mode is diffusive along the zone diagonals. There is also 
a sound-mode like mode at higher energy. Inset: dispersion of the 
collective modes along the TW(ir, tt) (solid symbols) and TX(0, tt) 
(open symbols) directions. The lower energy modes show a q 2 de- 
pendence in both directions except at (sufficiently) low q. The high 
energy mode is sound like with dispersion uj — vq. 



of 1*2/(1 + nq 2 ) and G2/(l + Kg 2 ), respectively. 

The order parameters associated with the parallel ( Ap) and 
diagonal (A^) nematic phases are given by 



Ap = F 2 ^2(cosk x - cosk y ) ^4.c k ^ , 

k 

Ad = G 2 J^(2 sin k x sin k y ) (c k c k \ . (2) 



Solving the above mean-field equation, Eq. (|2), one can ob- 
tain the behavior of nematic order parameters as one varies 
the ratio of F 2 / G 2 The transition to the diagonal nematic or- 
dered state from isotropic liquid phase occurs above a critical 
value of interaction G 2 , and it is second order as a function of 
chemical potential. The competition between the diagonal and 
the parallel nematics leads to the suppression of both phases, 
while they coexist in a finite window of chemical potential. 1 
The longer range hopping, such as t', is ignored in this study, 
since qualitative features are not affected by a finite t'J^Q- 



To study the behavior of non-Fermi liquid near the QCP, 
we set our parameter as follows: F 2 No = 0.1 and G 2 Nq = 
0.196, where N — l/(2tn 2 ) is the density of states at the 
Fermi surface. With these parameters, the QCP occurs near 
H c /(2t) = —0.2. While the parallel nematic also occurs at 
smaller value of [i within the above parameter space, below 
we will present the results near the /i c , since the non-Fermi 
liquid behavior is realized only near the QCP associated with 
the diagonal nematic phase. We first consider the dynamical 
effective interaction between the quasi-particles. 

Susceptibility in Random Phase Approximation (RPA) — 
The effective interaction mediated by the collective modes at 
RPA level is given by 



Xiv(q, v) 



G 2 (q) 



(3) 



l-G 2 (q)n° N (q^)' 

and the bare dynamical polarizability of the diagonal nematics 
is given by 



(27T) 
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C 2 (P) 2 , (4) 



where /(e) is the Fermi function. Near the nematic transition, 
the denominator in Eq. ([3]l becomes very small for q — > and 
v — ► 0, if v vanishes faster than q. This results in a singular 
behavior of the effective interaction in the forward scattering 
limit, and eventually non-Fermi liquid behavior. 

We carried out numerical integration of the above equa- 
tion. In Fig. [TJ the imaginary part of the effective interaction 
for the diagonal nematic is shown for small values of q for 
p,/2t = —0.22 slightly above the QCP. The diagonal nematic 
fluctuations has a propagating mode along the zone parallel 
directions and the mode becomes diffusive along the zone di- 
agonals (lmuj/uj is peaked at uj = 0). There exists also a 
sound-like mode at higher energy along the zone diagonals. 
These can be contrasted with the nematic Goldstone modes in 
the continuum case42i2i 

In the inset of Fig.Q] we plot the dispersion of these modes. 
For the highly damped mode along the diagonal direction, the 
dispersion is determined by the peak position of the imagi- 
nary part of the effective interaction, while the mode disper- 
sion along the zone parallel directions is given by the solution 
of 1 — G 2 (q)rio(q,cj(q)) = 0. The lower energy modes show 
a q 2 dependence in both directions, except at very low q val- 
ues with a finite but small gap of A/2t = 0.01. This mode 
becomes soft and the effective interaction becomes singular 
at the QCP. However, the numerical integration becomes non- 
trivial when q — > where v vanishes faster than q within the 
numerical computation. 

Quasiparticle Self-energy- To capture this singular behav- 
ior of the effective interaction, we expand the IT^ in the small 
q and small co/q limit as follows. 



IT 



bq 2 



ic- + 0[ — 

q \q z 



(5) 



where coefficients a, b and c can be determined numerically, 
and depend on the underlying band structure^i^ and tem- 
perature. The closeness to the QCP can thus be quantified 
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FIG. 2: (Color online) Calculated frequency dependence of the 
quasiparticle self-energy at the Fermi surface along (7r,7r) direction 
at different chemical potential, fj, — —0.3, —0.22 correspond to the 
symmetric phase, and /j, — —0.18 is inside the diagonal nematic 
phase, n = —0.2 is the closest point to the QCP and shows the 
strongest deviation from the Fermi liquid behavior at ui — > 0. Non- 
fermi liquid behavior is observed in all cases except at fi — —0.3, 
where Fermi liquid behavior is restored. The red (light) curves show 
the contribution from the regular part of IIjv to the self-energy. All 
quantities are in units of 2t. 

by quantity 6=1 — C?2(0)a. To understand how the de- 
cay rate for single-particle excitations is modified by the col- 
lective fluctuations, we calculate the imaginary part of the 
self-energy. At one-loop level, the imaginary part of the self- 
energy in the real frequency is given by 

ImE < k '">-/s/(S5&««k + f)x 

[b(v) + f(u + w)] Im X jv(q, v) Im G(k + q, uj + v) (6) 

In the lowest order approximation, we use the non-interacting 
Green's function and ImGo(k, ui) = ~tt5(lu — eu) in Eq. ||6). 
Figure |2] shows the imaginary parts of the self energies near 
the QCP both in the symmetric phase(^/2t = -0.3, -0.22) 
near the QCP (fi/2t = —0.2) and inside the diagonal ne- 
matic phase (fi/2t = —0.18) at the Fermi surface along the 
(tt, 7t) direction in the low T limit. At small ui, the self-energy 
shows clear deviation from the Fermi liquid behavior with the 
asymptotic behavior ImE ~ < 1. The scattering 

rate is the largest close to the QCP at \i]2t = —0.2, and 
self-energy shows a strong non-Fermi liquid behavior with 
a ~ 2 / 3 23,26 . This indicates that the quasiparticle excitations 
are not well defined approaching the Fermi surface, which 
rules out a Fermi-liquid description. This non-Fermi liquid 
behavior exists at all parts of the Fermi surface except at four 
points along the x and y axes, where the quasiparticle is well 
defined. In Fig. [2] we also compute the contribution from the 
regular part of 11^ to the self-energy, and we find the ampli- 
tude is always smaller than that of the singular part. This is 
due to the smaller density of states at the Fermi surface along 
the nodal directions. This should be contrasted to the paral- 
lel nematic order, where the density of the states at the Fermi 
surface is large due to the closeness of the Van Hove singular- 
ity near (tt, 0). The Fermi liquid behavior is restored far from 



FIG. 3: Calculated frequency dependence of the quasiparticle self- 
energy at the Fermi surface along (tt, tt) direction at different tem- 
peratures with fj, = —0.2. All quantities are in units of 2t. 



the QCP (/i = —0.3), where the singular and the regular parts 
overlap at low ui. 

In Fig. |3] we show the imaginary part of the self-energy at 
the Fermi surface in the nodal direction with fj,/2t = —0.2 
at different temperatures. This corresponds to the quantum 
critical regime above the QCP. 27 ' 28 Finite temperature effects 
enter through the finite temperature bose and fermi functions, 
as well as the temperature dependence of the expansion co- 
efficients in 11^. In addition to the singular contribution of 
IT^, the regular part is also included in the calculation. In 
the quantum critical regime at T > 0, the self-energy consists 
both contributions from the classical and the quantum fluc- 
tuations. The quantum part of the self-energy obeys (u>/T) 
scaling in the quantum critical regime, while the classical part 
dominates at uo — 0, and scales as Im E(kF,0) oc T£(T) oc 
\JT I InT. 23 ' 27 We observe this scaling of the classical fluctu- 
ations in the Im £(k.F, 0); however, our numerical precision 
can not separate these two contributions. 

Conclusion — A consequence of the nematic order is the 
distortion of the underlying Fermi surface due to the devel- 
opment of the nematic phase. It has been proposed that the 
Fermi surface distortion and critical Fermi surface fluctua- 
tions play an important role in strongly correlated electron 
systems. Recent experiments in Sr3Ru207 near the metam- 
agnetic QCP suggest that the formation of a spin-dependent 
nematic phase can explain the experimental observation of the 
strong anisotropy in the magnetoresistivities. 16,29 It has been 
also suggested that the in-plane anisotropy of spin dynamics 
in detwinned YBa2Cu30y_5 have indicated a possible exis- 
tence of two dimensional anisotropic liquid crystalline phase 
in high temperature cuprates i 14,15 ' 30 Close to the Pomeranchuk 
instability of the electronic system, the dynamics of the lat- 
tice distortion can be dramatically amplified. It has been ar- 
gued that electron-phonon interaction indeed plays an impor- 
tant role in the high-T c superconductors. 31 The anisotropic 
form factors of the B\ g buckling phonons has the form of 
(cosfcjc — cosfcy)^". Therefore, it is interesting to note that 
the effective interaction between electrons after integrating 
out the phonon mode has the form of d (k) in our Eq. ([TJ for 
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the small q limit. While the existence of a nematic phase still 
need further investigation, it is tempting to speculate that the 
non-Fermi liquid behaviors in these systems can be partially 
attributed to the closeness to the nematic-isotropic transition, 
where the Fermi surface becomes soft. 

In summary, we have studied the fluctuation effects of the 
diagonal electronic nematic order in the square lattice, using a 
model Hamiltonian with a quadrupole density-density interac- 
tion. We found that the collective mode associated with the di- 
agonal nematic fluctuations has a strong anisotropy with diffu- 
sive peaks around the Fermi surface except along the (0, ±7r) 
and (±7r, 0) directions. Close to the QCP, this mode becomes 
critical and strongly influences the decay rate of the quasi- 
particle; the imaginary part of the self-energy shows a non- 
Fermi liquid behavior. We found the decay rate at the Fermi 



surface is highly anisotropic and reaches a maximum at the 
QCP. It will be interesting to explore the effects of the collec- 
tive mode on charge and heat transports, Raman spectroscopy, 
and electron-phonon interaction near the transition. 
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